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Abstract. Turbulent accretion disks are the most probable origin of jets 
from young stellar objects. Using an extended version of ZEUS-3D code in 
the axisymmetry option we investigate the evolution of a disk wind into a 
collimated jet. Resistive jets propagate slower with increasing magnetic 
diffusivity, and the degree of collimation decreases. The resistive disk 
included in the simulations ejects persistent, well collimated jet. 



1. Magnetic jets 

Astrophysical jets are observed in objects with an accretion disk. The similarity 
of jets from the different sources implies that the basic mechanism for the jet 
formation should be the same. 

The mechanism which turns the inflowing, accreting matter into an disk 
outflow is still not understood. However, we know that for lifting the matter 
into the corona magnetic fields play a major role. In order to model the time- 
dependent evolution of jet formation we solve resistive MHD equations: 
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We perform two types of simulations: 1) The disk is defined as time- 
independent boundary condition (Ouyed &: Pudritz 1997; Fendt & Elstner 2000; 
Fendt & Cemeljic 2002). Here we investigate the acceleration and collimation of 
a resistive jet. 2) The evolution of the disk structure is included in the numerical 
simulation. 



2. Formation of resistive jets 

Within a subgrid of 280x80 grid points and physical grid of (ZxR)=(140x40)Ri 
we investigated the transition from accretion of matter in the resistive disk to 
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Figure 1. Global evolution of the jet on a grid of (ZxR)=(280x40)Ri 
with a resolution of 900 x 200 elements for different magnetic diffusivity. 
77=0,0.01,0.1 for top, middle and bottom panel, respectively. Colors 
(blue to yellow in decreasing manner) indicate density and lines denote 
linearly spaced poloidal field lines. The bow shock advances slower 
with increasing magnetic diffusivity. 
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Figure 2. Components of the velocities and time evolution of the 
mass flow ratio. Left panel: Cuts in Z direction at R=15 for the qua- 
sistationary time of simulations. Poloidal components are given in solid 
lines, toroidal in dashed lines (note the different left and right labels). 
Right panel: The radial/axial mass flow in the inner part of the jet, 
(ZxR)=(60x20)Ri. Parameter ?7=0, 0.1,0.5 in solid, dashed and dotted 
line, respectively. The degree of collimation decreases with increasing 
V- 

ejection in the jet for up to 4000 rotations of the inner disk radius Ri. The 
magnetic diffusivity, in combination with the magnetic and inertial forces and 
pressure and gravity, modifies the MHD structure of the jet. In our simulations 
the jet velocity inccreases with increasing magnetic diffusivity. The same holds 
for the degree of collimation. We defined the mass flux across the surfaces 
parallel to the accretion disk boundaries: 

Mz = / 2'KrpUzdr , = / 2'KrpUrdz . 
Jo Jo 

As a measure of the collimation we propose the ratio of the radial to the axial 
mass flow, where for this ratio >1 the flow becomes decollimated. We found 
that below the 'critical' r) (in our choice of parameters ij < 0.5), the outflow is 
collimated. 

The toroidal field may lead to (de-)accelerating Lorentz forces Fl y ~ j± x 
and (de-)collimating forces Fl,± ^ j\\ x -B (projections with respect to the 
poloidal magnetic field), see Figure 3. 

As a de-coUimation of the poloidal magnetic structure also implies a smaller 
launching angle for the sub-Alfvenic flow, the magneto-centrifugal acceleration 
mechanism may work more efficiently. 

3. Jet &Z. disk evolution 

The vertical disk equilibrium and the Keplerian rotation imposes the constraint 
for the radial density profile. The initial poloidal velocity profile is defined only 
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Figure 3. Lorentz forces in the jet for different magnetic diffusivity 
7] = 0, 0.1, 0.5 (solid, dashed and dotted hnes, respectively) along a 
flux surface leaving the box close to the (R=20,Z=60)-corner. For the 
perpendicular component the positive sign denotes the R-direction (de- 
coUimating force) . For the parallel component the positive sign denotes 
the Z-direction (accelerating force). 



by the accretion rate in the disk. The magnetic field is vertically uniform and 
the magnetic diffusivity is of a Gaussian profile. All the other boundaries are 
defined as open. 

Close to the origin, a 'sink' region is defined, ensuring the outflow to contain 
only the matter accreted from the disk, and not re-bounced from the symmetry 
axis of the jet (Casse Sz Keppens 2002) . 

Our preliminary results are as follows (Figure 5). The jet launching occurs 
only within R<5, indicating the existence of the conditions for the ejection only 
in the inner part of the jet (plasma /3 and angle of the magnetic field lines 
decrease for increasing radii), in agreement with Casse & Keppens (2002). The 
outflow is well collimated, with the velocity equal to the escape velocity of the 
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Figure 4. The density and the poloidal velocity. Cut in Z direction 
at R=15. 
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Figure 5. Evolution of the outflow from the accretion disk on a grid 
of (ZxR)=(60x40)Ri in a mesh of 125x80 elements for t=0,50,150 
rotations of the disk inner radii. Colors (blue to yellow in decreasing 
manner) indicate density and lines are linearly spaced poloidal field 
lines. 



central object. There exists strong indication for stationarity of the outflow after 
few tens of rotations. The jet radius is R~20. The disk piles-up for a factor two 
to three in comparison to the initial height, reaching new equilibrium. 
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